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Abstract 

The  effects  of  climate  change  are  expected  to  impact  military  installations 
in  the  Continental  United  States  (CONUS),  including  Army  installations 
that  have  large  land-based  range  areas  used  for  testing,  training,  or 
maneuvers.  Climate  change  has  the  potential  to  affect  several  management 
concerns  at  Army  installations.  Natural  areas  may  shift  on  installations 
and  change  the  costs  to  maintain  training  and  testing  areas.  Climate 
change  is  likely  to  increase  the  management  costs  for  Threatened  and 
Endangered  Species  (TES)  and  noxious  invasive  species  (NIS).  This 
document  describes  a  set  of  climate  change  data  gathered  to  support  a 
larger  project  undertaken  to  determine  the  thresholds  of  climatic 
characteristics  of  these  targeted  species.  This  work  describes  available 
climatic  data  that  will  show  what  the  major  Global  Climatic  Models 
(GCMs)  predict  about  those  changes,  using  Fort  Bragg,  NC  as  an  example. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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1  Introduction 

1.1  Background 

The  2007  Fourth  Assessment  Report  by  the  Intergovernmental  Panel  on 
Climate  Change  (IPCC  2007a)  states  that,  during  the  20th  century,  global 
surface  temperature  increased  0.6  ±  0.2  °C  (IPCC  2007b).  Much  of  the  ob¬ 
served  temperature  increase  since  the  middle  of  the  20th  century  has  been 
caused  by  increasing  concentrations  of  greenhouse  gases  (GHGs),  which 
result  from  human  activities  such  as  the  burning  of  fossil  fuel  and  deforesta¬ 
tion.  Global  dimming,  a  result  of  increasing  concentrations  of  atmospheric 
aerosols  that  block  sunlight  from  reaching  the  Earth’s  surface,  has  partially 
countered  the  effects  of  warming  induced  by  GHGs.  Climate  model  projec¬ 
tions  (IPCC  2007b)  indicate  that  global  surface  temperature  is  likely  to  rise 
between  1.1  and  6.4  °C  during  the  21st  century. 

The  President’s  Council  on  Environmental  Quality  (CEQ)  issued  draft 
guidance  (OSD  2010)  to  all  Federal  agencies  concerning  the  manner  in 
which  climate  change  should  be  included  in  the  evaluation  of  environmen¬ 
tal  effects  under  the  National  Environmental  Policy  Act  (NEPA).  Under 
NEPA,  Federal  agencies  are  required  to  evaluate  the  environmental  im¬ 
pacts  of  proposed  Federal  actions  and,  wherever  possible,  to  explore  a 
broad  range  of  options  for  minimizing  potentially  adverse  outcomes  and 
consequences  that  are  caused  -  wholly  or  in  part  -  by  those  actions.  This 
new  guidance  extends  the  issues  to  be  considered  to  include  GHG  emis¬ 
sions  and  climate  change,  and  how  agencies  should  address  the  interac¬ 
tions  between  their  proposed  actions  and  these  factors.  Specifically,  the 
guidance  states  that: 

With  regard  to  the  effects  of  climate  change  on  the  design  of  a  proposed 
action  and  alternatives,  Federal  agencies  must  ensure  the  scientific  and 
professional  integrity  of  their  assessment  of  the  ways  in  which  climate 
change  is  affecting  or  could  affect  environmental  effects  of  the  proposed 
action ... 


Climate  change  can  increase  the  vulnerability  of  a  resource,  ecosystem,  or 
human  community,  causing  a  proposed  action  to  result  in  consequences  that 
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are  more  damaging  than  prior  experience  with  environmental  impacts  anal¬ 
ysis  might  indicate ... 

Agencies  should  consider  the  specific  effects  of  the  proposed  action  (in¬ 
cluding  the  proposed  action’s  effect  on  the  vulnerability  of  affected  eco¬ 
systems),  the  nexus  of  those  effects  with  projected  climate  change  effects 
on  the  same  aspects  of  our  environment,  and  the  implications  for  the  en¬ 
vironment  to  adapt  to  the  projected  effects  of  climate  change  ... 


Where  agencies  consider  climate  change  modeling  to  be  applicable  to 
their  NEPA  analysis,  agencies  should  consider  the  uncertainties  associat¬ 
ed  with  long-term  projections  from  global  and  regional  climate  change 
models  ... 

As  with  other  agencies,  the  effects  of  climate  change  are  expected  to  impact 
Continental  United  States  (CONUS)  military  installations.  In  particular, 
Army  installations  have  large  land-based  range  areas  used  for  testing,  train¬ 
ing,  or  maneuvers.  Climate  change  has  the  potential  to  affect  several  man¬ 
agement  concerns  at  Army  installations.  Work  in  this  area  commonly  aimed 
at  the  management  of  Threatened  and  Endangered  Species  (TES)  and  the 
potential  appearance  and  increase  of  noxious  invasive  species  (NIS). 

1.2  Objective 

The  objective  of  this  effort  is  to  survey  available  Global  Climatic  Model 
(GCM)  data  and  access  the  set  of  “best  available  data”  intended  to  support 
the  larger  issue  of  determining  the  climatic  thresholds  for  various  plant 
and  animal  species. 

1.3  Approach 

Chapter  2  broadly  reviews  climate  change  research,  with  a  focus  on  the 
predicted  spatial  distribution  of  expected  changes. 

Chapter  3  describes  significant  datasets  and  the  procedures  found  in  a 
survey  of  climate  change  data.  GCM  outputs  of  temperature  and  precipita¬ 
tion  are  the  standard  metrics  likely  to  be  of  interest  to  TES  and  NIS  re¬ 
searchers,  both  areas  of  considerable  concern  to  Army  installation  land 
managers.  Other  researchers  (Busby  1991,  WorldClim  2012,  CCAFS  2011) 
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in  the  species  and  climatic  characteristics  fields  have  developed  a  set  of 
useful  derivative  metrics  that  were  found  to  be  of  great  use  in  this  project. 

Chapter  4  describes  in  detail  the  characteristics  of  the  best  data  found  and 
extracted  for  this  research,  and  provides  a  few  simple  illustrations  of  po¬ 
tential  applications.  Time  horizons  for  these  datasets  were  divided  into 
four  periods  spanning  from  1990  (the  base  year)  to  2085,  and  that  data 
used  in  illustrative  applications  is  specifically  for  Fort  Bragg,  NC. 

1.4  Scope 

This  investigation  reviews  the  available  literature  that  specifically  supports 
the  spatial  distribution  of  climatic  change  predictions.  This  work  made  no 
attempt  to  generate  new  predictions. 

1.5  Mode  of  technology  transfer 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW) 
at  URLs: 


http://www.cecer.armv.mil 

http://libweb.erdc.usace.armv.mil 
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2  Climate  Change  Modeling  Review 

2.1  General  background  to  climate  modeling 

Research  in  the  discipline  of  climate  change  dates  back  to  the  1960s.  the 
work  of  many  individuals  and  groups  to  objectively  understand  the  direc¬ 
tion  of  climate  change  has  generated  many  computer-based  models.  Some 
of  the  best  known  groups  include  National  Center  for  Atmospheric  Re¬ 
search  (NCAR,  in  Boulder,  CO),  the  Geophysical  Fluid  Dynamics  Labora¬ 
tory  (GFDL,  in  Princeton,  NJ),  the  Hadley  Centre  for  Climate  Prediction 
and  Research  (in  Exeter,  UK),  the  Max  Planck  Institute  for  Meteorology 
(in  Hamburg,  Germany),  and  the  Institut  Pierre-Simon  Laplace  (IPSL,  in 
Paris,  France). 

All  the  respected  models  generate  predictions  based  on  a  set  of  conven¬ 
tions  disseminated  through  the  IPCC.  Such  standardization  is  meant  to 
facilitate  comparison  between  models.  As  the  predictive  capabilities  of 
climatic  models  are  refined,  discrepancies  between  them  grow  less  signifi¬ 
cant.  However,  enough  variation  still  exists  that  critics  can  use  differences 
between  the  models  to  exaggerate  the  differences  within  climatic  research. 
To  minimize  such  confusion,  the  IPCC  acts  as  a  coordinating  organization 
and  its  reports  are  intended  to  reflect  the  scientific  consensus  of  the  ex¬ 
perts  in  the  field. 

2.2  Scenarios  upon  which  climate  modeling  efforts  are  based 

One  of  the  primary  responsibilities  of  the  IPCC  is  the  arrangement  of  a 
series  of  standard  future  scenarios  to  assist  with  coordination  and  compar¬ 
ison  between  model  inputs  and  their  results.  This  international  standard 
set  of  scenario  types  is  named  after  The  Special  Report  on  Emissions  Sce¬ 
narios  (SRES).  The  SRES  was  prepared  by  the  IPCC  for  the  Third  Assess¬ 
ment  Report  (TAR)  in  2001  on  future  emission  scenarios  to  be  used  to 
drive  GCMs  to  develop  climate  change  scenarios.  The  SRES  were  also  used 
for  the  Fourth  Assessment  Report  (AR4)  in  2007.  Table  1  lists  the  four 
scenario  families. 
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Table  1.  The  four  SRES  scenario  families  of  the  Fourth  Assessment  Report  with  associated 
projected  global  average  surface  temperature  increase  by  2099. 


Homogenous:  Global* 

Al 

Rapid  economic  growth 
(includes  groups:  AIT;  AIB;  AiFl) 

+1.4  -  6.4  °C 

BI 

Global  environmental  sustainability 

+1.1  -  2.9  °C 

Heterogeneous:  Regional  / 
Local 

A2 

Regionally  oriented 
economic  growth 

+2.0-5.4  °C 

B2 

Local  environmental  sustainability 

+1.4- 3.8  °C 

'Table  format  drawn  partially  from  IPCC  (2011). 

To  put  these  scenarios  in  a  slightly  different  light: 

Ai:  Maximum  energy  requirements  —  emissions  differentiated 

dependent  on  fuel  sources 
AiFl:  Fossil  intensive 

AiT :  Technology  development  of  non-fossil  sources 

AiB:  Balance  across  sources 

Bi:  Minimum  energy  requirements  and  emissions 

A2:  High  energy  requirements  —  emissions  less  than  AiFl 

B2:  Lower  energy  requirements  —  emissions  greater  than  Bi. 

2.3  The  major  climate  models 

Since  the  1990s,  the  international  climate  change  science  community  has 
participated  in  a  series  of  efforts  (often  called  “campaigns”)  to  carry  out  ma¬ 
jor,  mostly  coordinated  attempts  to  exercise  their  best  available  modeling 
capabilities  under  similar  sets  of  SRES  scenarios.  This  study  used  the  most 
recent  model  results,  i.e.,  the  AR4  (IPCC  2009).  Table  2  lists  the  major  par¬ 
ticipants  in  the  AR4  campaign  and  the  name  of  their  most  recent  models.  The 
next  major  coordinated  modeling  effort  will  be  the  IPCC  Fifth  Assessment 
Report  (AR5),  which  will  be  finalized  in  2014.  Modeling  efforts  for  AR5  have 
already  begun. 

Figure  1  shows  different  model  and  scenario  results  and  how  they  diverge 
over  time.  Near-term  results  can  be  considered  more  creditable  than  those 
that  approach  2099. 
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Table  2.  SRES  scenario  runs  for  AR4  (August  2006  data). 


Country 

Acronym 

Model  Name 

Beijing  Climate  Center 

China 

BCC 

CM1 

Bjerknes  Centre  for  Climate  Research 

Norway 

BCCR 

BCM2.0 

Canadian  Center  for  Climate  Modelling  and 
Analysis 

Canada 

CCCma 

CGCM3  (T47  resolution) 
CGCM3  (T63  resolution) 

Centre  National  de  Recherches 
Meteorologiques 

France 

CNRM 

CM3 

Australia's  Commonwealth  Scientific  and 
Industrial  Research  Organization 

Australia 

CSIRO 

Mk3.0 

Max-Planck-lnstitut  for  Meteorology 

Germany 

MPI-M 

ECHAM5-0M 

Meteorological  Institute,  University  of  Bonn 

Germany 

MIUB 

ECHO-G 

Meteorological  Research  Institute  of  KMA 

Korea 

METRI 

Model  and  Data  Groupe  at  MPI-M 

Germany 

M&D 

Institute  of  Atmospheric  Physics 

China 

LASG 

FGOALS-gl.O 

Geophysical  Fluid  Dynamics  Laboratory 

USA 

GFDL 

CM2.0 

CM2.1 

Goddard  Institute  for  Space  Studies 

USA 

GISS 

AOM 

E-H 

E-R 

Institute  for  Numerical  Mathematics 

Russia 

INM 

CM3.0 

Institut  Pierre-Simon  Laplace 

France 

IPSL 

CM4 

National  Institute  for  Environmental  Studies 

Japan 

NIES 

MIR0C3.2  hires 

MIR0C3.2  medres 

Meteorological  Research  Institute 

Japan 

MRI 

CGCM2.3.2 

National  Centre  for  Atmospheric  Research 

USA 

NCAR 

PCM 

CCS  M3 

UK  Met  Office  (Hadley  Centre) 

UK 

UKMO 

HadCM3 

HadGEMl 

National  Institute  of  Geophysics  and 
Volcanology 

Italy 

INGV 

SXG  2005 

*AR4  table  from  IPCC  (2010). 
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Figure  1.  Variation  between  GCMs  showing  ranges  within  scenarios. 
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2.4  Which  predictive  models  to  choose  for  this  work? 


To  determine  the  “best”  of  all  the  available  models  for  future  predictions  of 
different  variables,  one  objective  criterion  to  consider  is,  which  models 
have  demonstrated  the  closest  validation  with  the  variables  of  interest  — 
for  example,  precipitation  over  CONUS.  Those  models  that  have  had  the 
greatest  number  of  validation  studies  and  those  with  the  longest-period  of 
development  (l  to  2  decades)  include: 

1.  CM2.i  (GFDL  model  —  NOAA  Princeton) 

2.  E-H  and  E-R  (NASA  GISS) 

3.  HadGEMi  (Hadley  UKMO) 

4.  CGCM3  (Canadian  (CCCma)  model) 

5.  CCSM3(NCAR  Boulder). 


To  their  credit,  these  top  models  have  the  longest  history  and  largest  num¬ 
ber  of  peer  reviewed  publications.  Thus,  this  research  attempted  to  include 
as  many  of  these  models  as  possible,  and  to  exclude  models  with  shorter 
lifetimes  of  development,  and  fewer  person-hours  involved  in  their  valida¬ 
tion  and  calibration. 
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2.5  Problems  with  the  output  of  the  climate  models  and  the  desired 
solution 

The  characteristics  of  basic  climatic  models  presented  several  major  prob¬ 
lems  due  to.  First,  temperature  and  precipitation  (and  possibly  humidity) 
are  the  major  outputs  that  would  be  useful  for  the  characterization  of  spe¬ 
cies  changes.  Although  other  outputs  are  available,  they  do  not  relate  well 
to  issues  that  have  an  effect  on  the  life  habits  of  different  species.  Fortu¬ 
nately,  other  studies  (WorldClim  2012,  CCAFS  2011)  have  established  that 
it  is  possible  to  extract  data  that  are  highly  important  to  species  viability 
from  this  simple  temperature  and  precipitation  output.  Chapter  3  deals 
with  this  in  detail. 

The  second  issue  is  more  problematic.  GCMs  geographically  referenced 
data  are  gross  in  size.  In  fact,  most  GCMs  output  their  results  in  a  grid 
format  that  is  roughly  3x3  degrees  (-330  km  at  30  degrees  north)  in  size 
(see  Figure  2).  It  became  apparent  that  using  such  generalized  data  would 
result  in  less  than  satisfactory  results.  Fortunately  others  in  the  field  have 
agreed  that  this  was  an  issue  (Climate  Central  et  al.  2012,  CCAFS  2011) 
and  have  carried  out  “downscaling”  on  the  GCM  data.  “Downscaling,”  or 
refining  the  climate  model  results  to  specific  regions  involves  a  focus  on 
more  local  concerns  such  as  topography,  surface  winds,  evaporation,  and 
local  precipitation.  This  is  done  through  the  application  of  dynamic  and 
statistical  processes.  Statistical  downscaling  is  useful  for  adjusting  GCM 
scenarios  to  the  local  climate  observations  when  and  where  local  observa¬ 
tions  are  available  (Climate  Central  et  al.  2012).  Future  climate  scenarios 
has  been  downscaled  to  a  resolution  of  1/8  degree  (about  13  km,  see  Figure 
3)  using  statistical  approaches,  mostly  for  average  temperatures  and  pre¬ 
cipitation. 

Another  problem  was  the  nature  of  the  GCM  data.  At  the  project’s  outset, 
the  intention  was  to  begin  with  downscaled  data  acquired  for  a  previous 
effort  (Lozar  2011).  Figure  4  shows  a  range  of  temperature  predictions  for 
Fort  Bragg,  plotted  from  the  previously  compiled  data.  There  appear  to  be 
wide  temperature  swings  among  the  models  and  even  within  a  single  mod¬ 
el  in  the  first  few  decades  of  the  21st  century.  Later  (2050-2099),  the  varia¬ 
bility  seems  to  decrease,  but  this  could  easily  be  the  result  of  fewer  data 
points.  Both  the  A2  (red)  and  Bi  (green)  scenarios  seem  to  overlap  a  great 
deal,  such  that  no  pattern  appears  evident. 
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Figure  2.  The  coverage  of  a  single  temperature  data  point  from  Canadian  CGCM3 

model  (shown  in  green). 


Figure  3.  The  square  around  Fort  Bragg  is  3  degrees  on  edge-  one  GCM  pixel; 
much  greater  detail  is  made  available  by  downscaling  on  the  same  model  for  the 
year  2000  for  monthly  precipitation. 
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Figure  4.  Fort  Bragg  predicted  temperature  changes  for  seven  models  and  two  scenarios. 


After  reviewing  this  variable  GCM  data,  it  became  apparent  that  it  would 
be  difficult  to  justify  sensible  trends  based  on  data  that  varied  so  widely 
within  itself.  Discussions  with  John  W.  Weatherly  (ERDC-CRREL)  clari¬ 
fied  that  the  reason  behind  the  variation  was  that  the  data  were  weather 
predictions  rather  than  climate  predictions,  i.e.,  GCMs  are  designed  to 
mimic  weather  patterns  with  built-in  random  variation.  The  implication  of 
this  randomness  became  almost  immediately  apparent. 

To  extract  the  climatic  information  required  for  this  project  would  require 
averaging  the  data  over  a  20-  to  30-year  time  horizon,  using  the  midpoint 
date  as  the  climatic  data  point.  Such  an  effort  would  require  a  great  deal 
more  data  than  was  then  available,  plus  large  amounts  of  computer  calcu¬ 
lations  to  generate  the  desired  midpoint  data.  This  would  require  the  use 
of  several  years  of  data  on  either  temporal  side  of  the  target  year,  which 
would  be  averaged  to  develop  a  climatic  representation. 

Unfortunately,  such  averaging  removes  the  extreme  minimums  and  max- 
imums  that  have  the  greatest  effect  on  plant  and  animal  communities.  For 
non-averaged  data,  it  would  be  possible  to  use  the  downscaled  data  from 
Figure  4  that  previously  contained  too  much  randomness.  However,  ex¬ 
treme  weather  events  in  the  downscaled  data  are  not  predictive.  It  would 
be  possible  to  say  that  a  particular  maximum  temperature  would  be  N  °C 
and  that  in  too  years  it  would  occur  with  Y  frequency,  but  it  would  not  be 
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possible  to  say  accurately  when  those  critical  events  might  occur.  This 
seems  to  be  an  acceptable  tradeoff  for  those  who  model  thresholds,  but  it 
also  implies  it  would  be  necessary  to  download  all  too  years  of  downscaled 
data  instead  of  simply  using  data  from  a  few  snapshot  years.  Fortunately, 
others  (WorldClim  2012,  CCAFS  2011)  have  carried  out  the  required  data 
manipulation  to  generate  predictive  climatic  data  needed  for  this  work. 

Therefore,  the  ideal  result  of  this  effort  would  be  to  find  high-resolution 
downscaled  data  that: 

•  represents  more  than  just  temperature  and  precipitation 

•  has  been  averaged  over  a  multiple  decades  to  create  climatic  data 

•  would  provide  projections  in  the  relatively  near  term  (-2025s) 

•  include  at  least  five  established  GCMs  for  a  large  set  of  the  IPCC  sce¬ 
narios. 
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3  In  Search  of  the  Ideal  Dataset 

3.1  Preliminary  research 

The  revision  of  the  basic  dataset  was  begun  by: 

•  downloading  all  the  data  available  from  our  original  source 

•  reformatting  it  to  a  standard  GIS  form 

•  averaging  it  over  decades 

•  generating  results  similar  to  those  shown  in  Figure  4,  but  smoothed. 

At  the  beginning  of  this  process,  the  research  team  also  began  to  investi¬ 
gate  data  already  available  through  the  Internet,  particularly  that  available 
via  WorldClim.*  Climate  Change  Agriculture  and  Food  Security  (CCAFS 
2011)  uses  the  WorldClim  source  to  generate  19  corollary  biologically  re¬ 
lated  datasets,  and  these  additional  datasets: 

•  consecutive  months  -  the  maximum  number  of  consecutive  dry 
months  (<100  mm)  in  a  year 

•  precipitation  by  month 

•  temperature  maximum  by  month 

•  temperature  mean  by  month 

•  temperature  minimum  by  month. 

This  work  used  the  WorldClim  dataset,  which  represents  downscaled  data 
from  weather  stations  over  a  period  of  1950-2000, +  to  represent  “current” 
conditions. 

3.2  Cl  AT  data 

The  International  Centre  for  Tropical  Agriculture  (CIAT)  has  downscaled 
future  climate  projections  from  the  IPCC,  and  has  also  made  available  sets 
of  19  bioclimatic  concerns  (“bioclim”  data)  useful  in  characterizing  the 
biological  environment  and  the  predicted  GCM  changes.* 


*  http://www.worldclim.org/ 
t  Available  via  http://www.worldclim.org/current 
$  http://ccafs-climate.org/download  down.html 
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Figure  5.  CIAT  bioclimatic  data  content  illustration.  Where  Red  predominates  the  warmest 
month  is  significantly  warm.  Where  green  predominates,  the  coldest  month  is  very  cold. 
Where  Blue  predominates,  the  wettest  month  is  most  important  (e.g.,  the  rainforests  of  the 
Pacific  Northwest).  Similarly  darker  colors  suggest  the  area  is  lower  in  all  three  concerns 


bio_5_6_l  3CompositeBands3 
RGB 


■  Red:  Max  Temp  Warmest  Month 
I  I  Green:  Min  Temp  Coldest  Month 

■  Blue:  Prep  Wettest  Month 

These  are  comparable  at  30  arc-seconds  (~i  km)  resolution  to  the 
WorldClim  data  and  include  both  the  baseline  and  bioclimatic  data  layers. 
Since  this  work  is  interested  in  the  ecological  situation  at  a  specific  mili¬ 
tary  installation,  these  19  parameters  represent  many  of  the  concerns  that 
would  have  a  stressing  effect  on  the  flora  and  fauna  in  a  locality.  They  are 
derived  directly  from  the  base  temperature  and  precipitation  data  using 
ArcGIS  and  an  open  source  GIS  program,  DIVA-GIS  (Hijmans  2012).  The 
CIAT  web  site  clearly  documents  characteristics  of  the  data.  The  following 
paragraphs  describe  how  well  the  data  fulfills  the  stated  criteria  for  “ideal” 
high-resolution  downscaled  data. 

3.2.1  High-resolution  downscaled  data 

Although  CIAT  provides  several  different  resolution  datasets,  the  ~i-km2 
dataset  (their  highest  resolution)  was  selected.  These  data  are  downscaled 
to  a  similar  resolution  as  the  dataset  the  research  team  had  used  previous¬ 
ly  so  the  data  provide  some  variation  over  large  installations  and  within 
the  immediate  regions  of  any  Army  installation. 
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3.2.2  Averaged  over  a  few  decades  to  get  climatic  data 

The  CIAT  “Delta  method”  data  are  averaged  over  30  years  of  GCM  data 
with  the  average  centered  with  15  years  of  data  on  either  side  of  that  date. 
So,  the  data  meant  to  represent  the  “2025s”  is  an  average  of  GCM  output 
from  2010-2039.  Also,  the  WorldClim  site  provided  historical  data  that 
have  been  averaged  over  the  course  of  50  years  between  1950  and  2000. 

3.2.3  For  the  five  best  respected  GCMs  for  a  large  set  of  the  IPCC 
scenarios. 

CIAT  has  done  their  processing  on  many  GCMs,  including  the  entire  five 
target  GCMs  of  most  interest  to  this  project.  In  addition,  researchers  re¬ 
quested  a  sixth  GCM,  the  Australian  CSIRO  model,  to  increase  the  sample. 
Thus  the  adopted  GCMs  include: 

1 .  GFD  L  model  (N OAA  Princeton) 

6.  GISS  Model  e  (NASA  GISS) 

7.  UKMO  (UK  Hadley  Center) 

8.  CCCMA  (Canadian  model) 

9.  CCSM3(NCAR  Boulder) 

10.  CSIRO  (Australia). 

Note  that  three  scenarios  were  downloaded  for  each  GCM:  AiB,  A2,  and  Bi. 

3.2.4  From  which  could  be  extracted  some  reasonably  near-term  data 

With  the  baseline  data  from  the  20th  century,  the  time  horizon  and  inter¬ 
vals  for  the  data  used  are: 

•  late  20th  century,  or  “current” 

•  2010-2039  (henceforth  called  2025s) 

•  2040-2069  (henceforth  called  2055s) 

•  2070-2099  (henceforth  called  2085s). 

Thus  the  2025s  data  represent  a  time  frame  of  just  over  12  years  from  now, 
within  a  reasonable  planning  horizon  for  installations. 

3.2.5  Represented  bioclimatic  variables 

In  addition  to  the  basic  database,  WorldClim  also  provides  19  “bioclimatic” 
variables  that  have  been  derived  from  the  three  basic  climatic  variables 
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(minimum/maximum  temperature  and  precipitation),  and  that  are  found 
in  the  downscaled  GCM  data.  CIAT  has  calculated  the  same  bioclimatic 
layers  for  its  downscaled  future  projections,  and  has  provided  documenta¬ 
tion  for  deriving  the  bioclimatic  variables  from  any  basic  climate  change 
data  (Ramirez  and  Bueno-Cabrera  2009).  The  types  they  generated  are 
parallel  with  those  the  others  have  set  as  a  “Bioclimatic  Industry  Standard” 
(Phillips,  Anderson,  and  Schapire  2006). 

The  variables  included  in  the  CIAT  and  WorldClim  bioclimatic  data 
(Hijmans,  Cameron,  and  Parra  2012;  CCAFS  2011a)  are: 

•  Derived  from  maximum  and  minimum  temperature: 
o  BIOi  =  Annual  Mean  Temperature  (°C  times  10) 
o  BIO2  =  Mean  Diurnal  Range  (Mean  of  monthly 

(max  temp  -min  temp)) 

o  BIO3  =  Isothermality  (mean  diurnal  range/temperature 
annual  range) 

o  BIO4  =  Temperature  Seasonality  (standard  deviation  *100)  (°C 
times  10) 

o  BIO5  =  Max  Temperature  of  Warmest  Month 
o  BI06  =  Min  Temperature  of  Coldest  Month 
o  BIO7  =  Temperature  Annual  Range  (Bio5-Bio6) 
o  BI08  =  Mean  Temperature  of  Wettest  Quarter 
o  BIO9  =  Mean  Temperature  of  Driest  Quarter 
o  BIO10  =  Mean  Temperature  of  Warmest  Quarter 
o  BIO11  =  Mean  Temperature  of  Coldest  Quarter 

•  Derived  from  precipitation: 

o  BIO12  =  Annual  Precipitation  (in  millimeters) 
o  BIO13  =  Precipitation  of  Wettest  Month 
o  BIO14  =  Precipitation  of  Driest  Month 
o  BIO15  =  Precipitation  Seasonality  (Coefficient  of  Variation) 
o  BIO16  =  Precipitation  of  Wettest  Quarter 
o  BIO17  =  Precipitation  of  Driest  Quarter 
o  BIO18  =  Precipitation  of  Warmest  Quarter 
o  BIO19  =  Precipitation  of  Coldest  Quarter 

o  cons_mths  =  Consecutive  Months  -  the  maximum  number  of  con¬ 
secutive  dry  months  (<100  mm)  in  a  year. 
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3.3  Limitations  of  the  CIAT  data 

The  CIAT  data  reflect  the  standard  bioclimatic  datasets  being  generated. 
Still,  it  remains  to  be  determined— how  well  the  CIAT  data  fulfills  the 
needs  of  this  project,  i.e.,  the  management  of  TES  and  the  potential  ap¬ 
pearance  and  increase  of  NIS. 

3.3.1  Climate  vs.  weather  data 

A  primary  issue  in  TES  and  NIS  management  involves  changes  in  climatic 
parameters  that  pass  a  critical  threshold,  beyond  which  there  is  a  signifi¬ 
cant  difference  in  a  species’  viability  to  survive.  It  is  not  within  the  scope  of 
this  project  to  determine  the  thresholds  for  particular  species;  rather,  the 
objective  of  this  work  is  to  identity  the  climatic  data  that  supports  the 
identification  of  those  thresholds. 

Consultation  with  staff  working  on  related  research  (Hayden  2011)  re¬ 
vealed  that  thresholds  are  likely  to  be  found  at  the  extremes  of  climatic 
data  rather  than  at  the  norm  (norm  being  climatic  data,  extremes  to  be 
found  in  daily  weather  type  data).  Recall  that  the  problem  with  the  data 
from  a  previous  related  project  (Lozar  2011)  had  was  that  it  mimicked 
weather  data  with  a  built-in  random  component;  the  search  for  averaged- 
over-time  datasets  was  undertaken  to  overcome  this  problem.  It  may  seem 
reasonable  to  question  whether  it  is  possible  to  use  randomized  GCM  data 
for  climate  change  issues,  but  this  is  emphatically  not  the  case.  Because 
GCM  data  are  predictive,  any  occurrences  of  extreme  values  not  only  are 
fictitious,  but  random.  At  best,  one  could  say  that  an  extreme  value  occurs 
X  times  in  the  next  100  years.  However,  since  the  extreme  is  random,  one 
could  not  locate  that  event  in  time.  Therefore,  it  is  not  useful  to  know  the 
frequency  of  occurrence  that  could,  for  example,  occur  in  2012  or  2085. 

3.3.2  Forest  vegetation  simulator  categories 

Work  done  by  the  US  Forest  Service  at  the  Moscow  Forestry  Sciences  La¬ 
boratory  (USDA  2012)  created  data  that  would  be  useful  to  the  related 
threshold  project.  This  US  Forest  Service  date  were  similar  to  those  of  the 
CIAT  bioclimatic  derivations,  but  with  different  output  datasets  (created 
using  the  Forest  Vegetation  Simulator  [FVS]  modeling  system).  Unfortu¬ 
nately,  the  data  include  only  three  GCMs  for  the  years  2030,  2060,  and 
2090  and  only  extend  to  a  portion  of  the  Western  United  States  (and  thus 
does  not  include  Fort  Bragg).  Still,  it  is  instructive  to  compare  the  two  da- 
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tasets  to  see  how  well  the  FVS  (Table  3)  compares  with  CIAT  data  in  sup¬ 
porting  TES/NIS  threshold  work.  The  primary  issue  for  threshold  work  is 
that  both  datasets  use  averaged  GCM  data.  Thus,  neither  deals  with  ex¬ 
treme  occurrences.  In  this  critical  respect,  they  are  similar. 

Note  that  “degree-days”  unit  in  Table  3  is  used  in  estimating  the  demand 
for  energy  required  for  heating  or  cooling.  In  the  United  States,  the  typical 
standard  indoor  temperature  is  65  °F.  For  each  1  °F  decrease  or  increase 
from  this  standard  in  the  average  outside  temperature,  1  heating  or  cooling 
degree  day  is  recorded.  For  example,  an  average  outside  temperature  for  a 
day  of  60  °F,  records  as  5  heating  degree-days  (HDD);  if  it  were  70  °F,  it 
would  record  as  5  cooling  degree-days  (CDD). 

These  concerns  focus  on  issues  relevant  to  questions  of  stress  on  vegeta¬ 
tive  communities.  Since  these  data  are  available  for  only  Western  US  loca¬ 
tions  at  this  time,  data  was  requested  for  Fort  Irwin,  CA  (Appendix  A).  Of 
the  18  variables  listed,  only  12  (highlighted  in  grey  in  Table  3)  could  be 
identified  in  the  Fort  Irwin  output. 


Table  3.  The  FVS  set  of  concerns. 


1 

Mean  annual  temperature 

2 

Average  temperature  in  the  coldest  month 

3 

Minimum  temperature  in  the  coldest  month 

4 

Average  temperature  in  the  warmest  month 

5 

Maximum  temperature  in  the  warmest  month 

6 

Annual  precipitation 

7 

Growing  season  precipitation:  April  -  September 

8 

Summer-winter  temperature  differential  (Variable  4  -  Variable2) 

9 

Degree-days  >5  0  C 

10 

Degree-days  <0  °C 

11 

Minimum  temperature  of  degree-days  <0  °C 

12 

Julian  date  of  the  last  freezing  date  of  spring 

13 

Julian  date  of  the  first  freezing  date  of  autumn 

14 

Length  of  the  frost-free  period 

15 

Accumulated  degree-days  >5  °C  within  the  frost-free  period 

16 

Julian  date  when  the  sum  of  degree-days  >5  °C  reaches  100 

17 

Annual  dryness  index:  Ratio  of  Variable  9  to  Variable  6 

18 

Summer  dryness  index:  Ratio  of  Variable  15  to  Variable  7 
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3.3.3  FVS  data  categories  compared  to  the  CIAT  categories 

A  side-by-side  comparison  of  FVS  and  CIAT  categories  of  resolves  the  ques¬ 
tion  of  which  of  the  desired  FVS  data  types  can  be  obtained  from  the  CIAT 
bioclimatic  data,  and  which  sets  are  directly  comparable.  Table  4  lists  eight 
categories  that  are  exactly  the  same  as,  or  that  can  be  generated  from,  the 
bioclimatic  data.  Most  of  the  remaining  FVS  categories  can  be  generated 
from  the  CIAT  datasets  with  varying  degrees  of  estimated  difficulty  (also 
indicated).  Thus,  the  CIAT  data  well  represent  the  FVS  categories. 

For  the  purposes  of  finding  species’  bioclimatic  thresholds,  the  questions 
that  immediately  follow  from  examining  Table  4  are: 

•  How  important  is  each  of  the  listed  concerns  (i.e.,  do  we  really  need 
them  all?) 

•  Are  some  critical  categories  and  need  to  be  calculated  /  generated? 

•  Are  there  critical  others  not  on  our  list? 

In  some  cases,  the  calculation  effort  is  easy,  but  the  importance  of  the  re¬ 
sult  is  not.  Also,  it  is  desirable  to  cover  any  situation  where  the  importance 
is  great.  To  be  efficient  and  effective,  it  is  also  desirable  to  avoid  situations 
where  the  effort  is  low,  but  the  payback  is  also  low.  So  the  next  step  is  to 
prioritize  (and  possibly  augment)  the  FVS  list  so  that  it  helps  to  directs  the 
next  efforts.  A  prioritized  list  will  be  available  after  work  on  this  portion  of 
the  project  ends  (Hayden  2011). 

3.4  The  concept  of  spatial  thresholds 

In  the  classification  of  zones  (whether  ecological,  social,  cultural,  or  oth¬ 
erwise),  the  center  of  a  zone  represents  the  greatest  degree  of  hardiness. 
However,  it  is  at  the  edges  of  the  zones  that  the  most  diverse  and  dynamic 
activity  takes  place.  Most  often  these  “edges”  are  blurry  transitions  where 
the  characteristics  of  two  or  more  classifications  mix  and  interact.  An  im¬ 
portant  example  are  estuaries,  which  occur  when  fresh  water  streams  en¬ 
ter  bodies  of  salt  water.  These  transition  areas,  which  are  among  the  most 
biologically  diverse  places  in  the  world,  are  absolutely  critical  to  the  health 
of  an  ecosystem.  With  this  in  mind,  the  climate  change  data  were  exam¬ 
ined  to  identify  and  highlight  those  areas  in  which  change  occurred  rapidly 
between  two  specific  classifications  (thresholds). 
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Table  4.  FVS  category  compared  to  CIAT  Bioclimatic  category  indicating  the  estimated 
Difficulty  to  Calculate  (5=Difficult,  the  same  data  show  the  Difficulty  =  0). 


FVS# 

FVS  Name 

CIAT# 

CIAT  Name,  or  Can  It  Be  Generated? 

Ease  to 
Generate 

1 

Mean  annual  temperature 

BI01 

Annual  Mean  Temperature 

0 

2 

Average  temperature  in  the 
coldest  month 

BI011 

Mean  Temperature  of  Coldest  Quarter 
(or  can  be  derived  from  monthly  data) 

3 

3 

Minimum  temperature  in 
the  coldest  month 

BI06 

Min  Temperature  of  Coldest  Month 

0 

4 

Average  temperature  in  the 
warmest  month 

BI010 

Mean  Temperature  of  Warmest 

Quarter  (or  can  be  derived  from 
monthly  data) 

0 

5 

Maximum  temperature  in 
the  warmest  month 

BI05 

Max  Temperature  of  Warmest  Month 

0 

6 

Annual  precipitation 

BI012 

Annual  Precipitation 

0 

7 

Growing  season 
precipitation:  April  - 
September 

Generate  (Sum  of  Precip  from  April- 
Sept) 

1 

8 

Summer-winter 
temperature  differential 
(Variable  4  -  Variable2) 

BI07 

Temperature  Annual  Range  (P5-P6) 

0 

9 

Degree-days  >5  °C 

Roughly  generate  from  monthly  values 
>5  °C  by  Month 

4 

10 

Degree-days  <0  °C 

Roughly  generate  from  monthly  values 
<0  °C  by  Month 

4 

11 

Minimum  temperature  of 
degree-days  <0  °C  (coldest 
temperature) 

Generated  by  finding  coldest 
temperature  by  month 

2 

12 

Julian  Date  of  the  last 
freezing  date  of  Spring 

Generated  by  taking  the  monthly 
trend  slopes’  intersection  with  0  °C 

5 

13 

Julian  Date  of  the  first 
freezing  date  of  Autumn 

Generated  by  taking  the  monthly 
trend  slopes’  intersection  with  0  °C 

5 

14 

Length  of  the  frost-free 
period 

Generate  from  Min  Temperature  by 
Month  or  FVS#13-12 

1 

15 

Accumulated  degree-days 
>5  °C  within  the  frost-free 
period 

Generate  from  selected  monthly 
temps  in  the  period  FVS#12-13 

3 

16 

Julian  Date  when  the  sum 
of  degree-days  >5  °C 
reaches  100 

Many  unclear  steps 

5 

17 

Annual  dryness  index: 

Ratio  of  Variable  9  to 
Variable  6 

Easy  to  calculate 

1 

18 

Summer  dryness  index: 

Ratio  of  Variable  15  to 
Variable  7 

Easy  to  calculate 

1 

ERDC/CERL  TR-13-23 


20 


In  one  of  the  derivative  bioclimatic  layers  (Bio  9  -  Mean  Temperature  of 
Driest  Quarter),  the  thresholds  between  classification  are  clearly  delineat¬ 
ed.  Figure  5  shows  the  temporal  and  spatial  migration  of  thresholds  be¬ 
tween  classifications:  cool  colors  (blues)  indicate  areas  of  stark  thresholds 
that  existed  in  the  second  half  of  the  20th  century.  Warm  colors  (reds,  yel¬ 
lows)  indicate  the  location  of  intense  thresholds  in  the  given  projected 
time  period  (the  2020s,  2050s,  and  2080s).  Green  areas  indicate  areas  of 
classification  stability. 

In  other  layers,  such  as  general  temperature  and  precipitation,  spatial 
shifts  occur  in  a  smooth  gradient.  Significant  thresholds  appear  only  in 
mountainous  areas  (i.e.,  where  the  threshold  lines  are  close  together  in 
Figure  6).  The  location  and  shift  of  spatial  thresholds  can  easily  be  calcu¬ 
lated  for  minimum/maximum/mean  precipitation,  temperature  and  any 
of  the  19  derived  bioclimatic  variables  described  in  Section  3.2.5  . 


Figure  6.  Bio  9  threshold  migration.  Cool  colors  (blues)  indicate  sharp  thresholds 
in  the  20th  century;  warm  colors  (reds)  indicate  intense  thresholds  in  each  of  the 
three  projected  future  periods. 
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4  Dataset  Characterizations 

4.1  Fort  Bragg  data 

The  centroid  climatic  parameter  from  the  CIAT  dataset  was  extracted  for 
all  the  data  types  for  Fort  Bragg,  NC.  (Appendix  B  to  this  report  includes 
the  complete  dataset.) 

4.2  What  the  CIAT  Bragg  data  include 

This  work  selected  the  following  parameters  for  the  evaluations  dealing 
with  installations: 

•  Limited  to  the  installation  Fort  Bragg  (for  point  centroid  not  including 
Camp  MacKall  Military  Reservation).  The  same  data  can  be  accessed 
for  anywhere  in  the  world. 

•  Limited  to  four  time  snapshots: 
o  Baseline  (1950-2000) 

o  2010-2039  (the  2020s) 
o  2040-2069  (the  2050s) 
o  2070-2099  (the  2080s) 

•  Using  results  from  six  GCMs: 

o  GFDL  model  (NOAA  Princeton):  gfdl_cm2_i 
o  GISS  Model  e  (NASA  GISS):  giss_model_er 
o  HadGEMi(UKMO):  ukmo_hadcm3 
o  CCCMA  (Canadian  model):  cccma_cgcm3_i_t47 
o  CCSM3  (NCAR  Boulder):  ncar_ccsm3_o 
o  CSIRO  (Australia):  csiro_mk3_5 

•  Limited  to  three  IPCC  scenarios: 

o  Scenario  Alb:  Balance  across  sources 
o  Scenario  A2:  Regional  Rapid  Economic  Growth 
o  Scenario  Bi:  Global  Environmental  Sustainability 

•  Using  68  BioClimatic  predictions: 

o  BIOi  =  Annual  Mean  Temperature  (deg  CX10) 
o  BIO2  =  Mean  Diurnal  Range  (Mean  of  monthly  (max  temp  -min 
temp)) 

o  BIO3  =  Isothermality  (P2/P7)  (*  100)  (Mean  diurnal 
range/temperature  annual  range) 
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o  BIO4  =  Temperature  Seasonality  (standard  deviation  *100)  (deg 
CX10) 

o  BIO5  =  Max  Temperature  of  Warmest  Month 
o  BI06  =  Min  Temperature  of  Coldest  Month 
o  BIO7  =  Temperature  Annual  Range  (P5-P6) 
o  BI08  =  Mean  Temperature  of  Wettest  Quarter 
o  BIO9  =  Mean  Temperature  of  Driest  Quarter 
o  BIO10  =  Mean  Temperature  of  Warmest  Quarter 
o  BIO11  =  Mean  Temperature  of  Coldest  Quarter 
o  BIO12  =  Annual  Precipitation  (mm) 
o  BIO13  =  Precipitation  of  Wettest  Month 
o  BIO14  =  Precipitation  of  Driest  Month 
o  BIO15  =  Precipitation  Seasonality  (Coefficient  of  Variation) 
o  BIO16  =  Precipitation  of  Wettest  Quarter 
o  BIO17  =  Precipitation  of  Driest  Quarter 
o  BIO18  =  Precipitation  of  Warmest  Quarter 
o  BIO19  =  Precipitation  of  Coldest  Quarter 

o  Consecutive  Months  -  the  maximum  number  of  consecutive  dry 
months  (<100  mm)  in  a  year, 
o  Precipitation  by  Month 

o  Temperature  Maximum  by  Month  (12  predictions) 
o  Temperature  Mean  by  Month  (12  predictions) 
o  Temperature  Minimum  by  Month  (12  predictions) 

4.3  Size  of  dataset  for  a  location 

For  the  installation  as  point  data,  5544  items/point  are  derived  from: 

•  68  climate  variables  (minimum/maximum/mean  temperature  and 

precipitation  x  12  months,  19  bioclimatic  variables,  one  consecutive  dry 
months  calculation)  for: 

o  six  GCMs  (GFDL,  GISS,  UKMO,  CCCMA,  CCSM3,  CSIRO),  for: 

*  three  scenarios  (SRES  AiB,  SRES  A2,  SRES  Bi),  and 

*  three  timeframes  (2010-2039  (the  “2025s”),  2040-2069  (the 
“2055s”),  and  2070-2099  (the  “2085s”). 

There  also  exist  in  the  data  the  same  68  climate  variables  for  the  second 
half  of  the  20th  century  (1950-2000)  downscaled  from  climate  station  data 
by  WorldClim  and  centered  on  1990.  (Appendix  B  includes  all  the  Fort 
Bragg  data.) 
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4.4  Examples  of  how  the  data  can  be  used. 

Figure  7  shows  how  the  change  in  the  quality  of  the  data  source  has  im¬ 
proved  the  data  quality  over  that  shown  in  Figure  4  (p  10).  Although  varia¬ 
tion  still  exists,  it  is  evident  that  the  data  in  Figure  7  are  much  more  direc¬ 
tional,  less  random,  and  more  reliable  for  their  intended  purpose  than  the 
previous  values. 

Figure  8  shows  the  most  basic  application  of  the  climatic  data  at  Fort 
Bragg,  NC.  The  left  graph  shows  all  the  GCM  predictions  for  the  “2025” 
time  frame.  Each  “grouping”  shows  the  effect  of  the  different  scenarios 
(AiB,  A2,  and  Bi).  Each  model  is  represented  by  a  differently  colored 
line.  The  right  graph  shows  the  same  data  for  the  next  period,  the  “2055” 
time  frame.  (Figure  8  does  not  show  changes  between  time  frames.) 

Applying  basic  statistical  manipulations  to  the  data  in  Figure  8  yields  the 
simpler  averaged  summary  (Figure  9).  This  average  of  all  the  models  indi¬ 
cates  that  the  temperature  will  steadily  rise  (black  line)  by  nearly  3  °C. 
Significantly,  in  the  near  term,  it  will  rise  by  a  full  degree  by  the  year  2025. 
The  standard  deviation  lines  show  that  there  is  good  agreement  among  all 
the  models  for  all  scenarios  in  the  near  term  (to  2025)  and  that  they  widen 
as  time  goes  on.  The  greatest  variation  between  the  maximum  and  mini¬ 
mum  predictions  is  about  3.5  °C.  All  the  lines  on  this  graph  indicate  a 
temperature  increase  not  less  than  1.5  °C  (current  to  2085  minimum  val¬ 
ue).  This  chart  reflects  the  expected  interpretation  of  GCM  climatic  data. 

One  of  the  CIAT  datasets  indicated  consecutive  dry  months  -  the  maximum 
number  of  months  during  the  year  with  less  than  100  mm  of  rainfall.  As  an¬ 
other  simple  example,  one  may  ask  whether  the  number  of  consecutive  dry 
months  between  the  baseline  and  the  2025s  will  change.  The  answer  to  this 
question  is  a  simple  matter  of  subtracting  the  2025s  data  (scenario  AiB  of  the 
Canadian  GCM)  from  the  baseline  data.  Figure  10  shows  the  projected 
changes  in  consecutive  dry  months  from  the  baseline.  In  Figure  10,  the  light¬ 
er  color  shows  greater  change.  White  represents  an  increase  of  4  consecutive 
dry  months  (Camp  MacKall).  Most  of  Fort  Bragg  will  experience  an  increase 
of  2  consecutive  dry  months  (mid-dark  grey).  It  is  likely  that  this  near-term 
change  will  be  of  interest  to  flora  and  fauna  threshold  research. 
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Figure  7.  CIAT  Fort  Bragg  temperature  (xlOO)  predictions;  cf.  Fig.  4. 


Figure  8.  CIAT  data  displayed.  On  the  left,  Fort  Bragg  temperatures  for  “2025s”;  on  the  right, 

compared  to  similar  2055s  data. 
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Figure  9.  Bio  1  —  Annual  mean  temperature  (°C  xlO):  average,  standard  deviation  (SD), 
maximum,  minimum,  all  GCMs  and  scenarios. 


Figure  10.  Increased  number  of  consecutive  dry  months  from  1950-2000  to  the  2020s  at 

Fort  Bragg. 
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5  Summary  and  Recommendations 

5.1  Summary 

Based  on  the  stated  needs  of  the  parent  project,  Prediction  and  Adaptation  of 
Military  Natural  Infrastructure  in  Response  to  Climate  Change,  this  work 
has  identified  the  most  extensive  bioclimatic  dataset  available  that  has  the 
potential  to  provide  the  data  needed  to  determine  how  climatic  change  might 
affect  thresholds  for  the  survival  and  migration  of  species  of  interest  to  mili¬ 
tary  land  managers.  These  data  categories,  described  in  Section  4.2  (p  21), 
were  generated  by  the  International  Center  for  Tropical  Agriculture,  CIAT, 
Cali,  Colombia,  and  are  now  available  to  support  the  parent  project. 

5.2  Recommendations 

This  work  has  provided  a  solid  foundation  for  the  parent  project  to  deter¬ 
mine  how  climatic  change  might  affect  thresholds  for  the  survival  and  mi¬ 
gration  of  species  of  interest  to  military  land  managers.  It  is  recommended 
that  future  work  continue  to: 

•  generate  data  for  specific  military  installations 

•  generate  spatial  data  over  a  region  centered  on  an  installation 

•  generate  new  bioclimatic  data  via  calculations  from  the  existing  set  of  data 

•  generate  spatially  related  threshold  maps  over  time  as  the  climatic  pa¬ 
rameters  shift  due  to  climatic  change. 

•  generate  the  likely  migration  of  focus  species  as  the  sister  work  on  this 
project  matures  in  determining  thresholds  per  species,  including  the: 
o  migration  of  focus  species  off  an  installation 

o  migration  of  focus  (likely  invasive  species)  onto  an  installation 

•  address  broader  concerns  specifically  related  to  military  specific  issues, 
such  as: 

o  viability  of  current  training  lands  over  time 
o  problems  of  increasing  erosion  on  military  lands 
o  identification  of  new  military  lands  in  the  CONUS  similar  to  focus 
lands  elsewhere  in  the  world 

o  long-term  viability  of  Military  missions  at  different  installations  as 
climatic  factors  change. 
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Acronyms  and  Abbreviations 


Term  Definition 

BCC  Beijing  Climate  Center 

BCCR  Bjerknes  Centre  for  Climate  Research 

CCAFS  Climate  Change  Agriculture  and  Food  Security 

CCCma  Canadian  Center  for  Climate  Modelling  and  Analysis 

CEERD  US  Army  Corps  of  Engineers,  Engineer  Research  and  Development  Center 

CEQ  Council  on  Environmental  Quality 

CERL  Construction  Engineering  Research  Laboratory 

CIAT  International  Center  for  Tropical  Agriculture  (CIAT,  Centro  Internacional  de 

Agricultura  Tropical,  Cali,  Colombia) 

CNRM  Centre  National  de  Recherches  Meteorologiques 

CONUS  Continental  United  States 

CSIRO  Australia's  Commonwealth  Scientific  and  Industrial  Research  Organization 

ERDC  Engineer  Research  and  Development  Center 

FVS  Forest  Vegetation  Simulator 

GCM  Global  Climate  Model 

GFDL  Geophysical  Fluid  Dynamics  Laboratory 

GFIG  Greenhouse  Gas 

GISS  Goddard  Institute  for  Space  Studies 

INGV  National  Institute  of  Geophysics  and  Volcanology 

INM  Institute  for  Numerical  Mathematics 

IPCC  Intergovernmental  Panel  on  Climate  Change 

IPSL  Institut  Pierre-Simon  Laplace 

LASG  Institute  of  Atmospheric  Physics 

M&D  Model  and  Data  Groupe  at  MPI-M 

METRI  Meteorological  Research  Institute  of  KMA 

MIUB  Meteorological  Institute,  University  of  Bonn 

MPI-M  Max-Planck-lnstitut  for  Meteorology 

MRI  Meteorological  Research  Institute 

NCAR  National  Center  for  Atmospheric  Research 

NEPA  National  Environmental  Policy  Act 

NIES  National  Institute  for  Environmental  Studies 

NIS  Noxious  Invasive  Species 

OSD  Office  of  the  Secretary  of  Defense 

SD  Standard  Deviation 

SF  standard  form 

SRES  The  Special  Report  on  Emissions  Scenarios 
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Term 

Definition 

TAR 

Third  Assessment  Report 

TES 

Threatened  and  Endangered  Species 

TR 

Technical  Report 

UK 

United  Kingdom 

UKMO 

UK  Met  Office  [Hadley  Centre] 

URL 

Universal  Resource  Locator 

USA 

United  States  of  America 

WWW 

World  Wide  Web 
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Appendix  A:  FVS  Climate 
Attributes  for  Fort  Irwin,  CA 


Table  Al.  FVS  climate  attributes  for  Fort  Irwin,  CA. 


Mean 

Annual 

Temp 

Mean 

Annual 

Precip 

Growing 

Season 

Precip 

Min  Temp 
Coldest 

Month 

Min  Temp  of 
degree-days 
<0  °C 

Max  Temp 
Warmest 

Month 

Summer-Winter 
Temp  Differential 

Frost-Free 

Period 

Degree- 

Days 

>5"C 

Accum 
Degree-Days 
>5  °C  within 
ffp 

Julian  Date 

When  the 

Sum  of 
Degree-Days 
>5  °C 

Reaches  100 

Degree-Days 

<0*C 

Scenario 

Year 

mat 

map 

gsp 

mtcm 

mmin 

mtwm 

mmax 

sday 

ffp 

dd5 

gsdd5 

dlOO 

ddO 

pSite 

CGCM3_A1B 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

CGCM3_A1B 

2030 

23.9 

87 

27 

11.6 

3.9 

35.9 

44.5 

12 

346 

6905 

6767 

13 

0 

0 

CGCM3_A1B 

2060 

24.9 

76 

25 

12.1 

4.5 

37 

45.5 

10 

357 

7270 

7171 

13 

0 

0 

CGCM3_A1B 

2090 

25.8 

73 

19 

13.2 

5.2 

38.4 

47 

12 

351 

7580 

7431 

11 

0 

0 

CGCM3_A2 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

CGCM3_A2 

2030 

23.6 

82 

28 

10.9 

3.5 

36.4 

45.2 

10 

341 

6801 

6673 

15 

0 

0 

CGCM3_A2 

2060 

25.3 

77 

34 

13 

4.9 

37.5 

46.1 

9 

357 

7378 

7274 

12 

0 

0 

CGCM3_A2 

2090 

27.1 

48 

18 

14.1 

5.8 

38.8 

47.3 

16 

346 

8037 

7845 

11 

0 

0 

CGCM3_B1 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

CGCM3_B1 

2030 

24 

76 

23 

11.4 

3.5 

36.9 

45.8 

23 

328 

6952 

6688 

12 

0 

0 

CGCM3_B1 

2060 

24.4 

73 

20 

12.2 

4.5 

36.8 

45.5 

17 

343 

7082 

6896 

13 

0 

0 

CGCM3_B1 

2090 

24.7 

77 

23 

12.5 

4.6 

37.4 

45.9 

9 

357 

7197 

7097 

13 

0 

0 

GFDLCM21_A2 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

GFDLCM21_A2 

2030 

23.9 

76 

23 

10.4 

2.6 

37 

45.7 

20 

329 

6903 

6698 

16 

0 

0 

GFDLCM21_A2 

2060 

25 

75 

19 

12.6 

4.6 

38.5 

47.1 

17 

342 

7293 

7087 

13 

0 

0 

GFDLCM21_A2 

2090 

27.3 

46 

17 

13.4 

5.1 

41.2 

49.9 

16 

348 

8115 

7919 

12 

0 

0 

GFDLCM21_B1 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

GFDLCM21_B1 

2030 

23.6 

68 

21 

11.3 

3.4 

36.4 

45.3 

30 

317 

6804 

6478 

15 

0 

0 

GFDLCM21_B1 

2060 

24.4 

64 

20 

12.1 

4.1 

37.1 

45.9 

23 

334 

7093 

6828 

13 

0 

0 

GFDLCM21_B1 

2090 

24.7 

75 

25 

11.3 

3.3 

38.1 

46.9 

21 

337 

7190 

6963 

13 

0 

0 

HADCM3_A2 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

HADCM3_A2 

2030 

23.6 

86 

30 

10.8 

2.7 

36.5 

45 

19 

334 

6792 

6598 

15 

0 

0 
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Mean 

Annual 

Temp 

Mean 

Annual 

Precip 

Growing 

Season 

Precip 

Min  Temp 
Coldest 

Month 

Min  Temp  of 
degree-days 
<0  °C 

Max  Temp 
Warmest 

Month 

Summer-Winter 
Temp  Differential 

Frost-Free 

Period 

Degree- 

Days 

>5"C 

Accum 
Degree-Days 
>5  °C  within 
ffp 

Julian  Date 

When  the 

Sum  of 
Degree-Days 
>5  °C 

Reaches  100 

Degree-Days 

<0*C 

Scenario 

Year 

mat 

map 

gsp 

mtcm 

mmin 

mtwm 

mmax 

sday 

ffp 

dd5 

gsdd5 

dlOO 

ddO 

pSite 

HADCM3_A2 

2060 

25.1 

133 

73 

12.1 

3.7 

38.9 

47.1 

21 

331 

7327 

7074 

14 

0 

0 

HADCM3_A2 

2090 

27.2 

114 

65 

13.2 

4.8 

41.5 

49.5 

13 

344 

8084 

7917 

12 

0 

0 

HADCM3_B2 

1990 

22.1 

76 

26 

9.9 

1.9 

34.5 

43.3 

41 

294 

6234 

5858 

19 

0 

0 

HADCM3_B2 

2030 

23.9 

90 

40 

11.4 

2.9 

36.5 

45.3 

27 

324 

6903 

6644 

15 

0 

0 

HADCM3_B2 

2060 

24.9 

116 

64 

11.8 

3.5 

37.9 

46.7 

21 

329 

7277 

7023 

13 

0 

0 

HADCM3_B2 

2090 

25.4 

140 

75 

13 

4 

38.2 

47 

15 

343 

7455 

7268 

12 

0 

0 

Note  that  the  data  distributed  has  abbreviations  as  column  headers.  The  extended  description  is  a  guess  as  to  what  the  abbreviations  mean.  The  web  site  lists  18  variables  as  output,  but  only  14  columns  of  climatic  parameters  are  present  in 
the  output.  Other  columns  were  present  in  the  output  associated  with  different  tree  species;  they  are  not  included  in  this  table. 
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Appendix  B:  CIAT  Dataset  for  Fort  Bragg,  NC 


Table  Bl.  CIAT  dataset  for  Fort  Bragg,  NC. 


Fort  Bragg  Centroid 

Theme 

Month 

Year 

Model 

Scenario 

1990 

2010-2039 

2040-2069 

2070-2099 

cccma_cgcm3_l_t47 

A1B 

bio 

1 

157 

164 

180 

178 

2 

126 

100 

126 

97 

3 

38 

33 

38 

31 

4 

7464 

7652 

7793 

7806 

5 

315 

307 

341 

321 

6 

-10 

7 

14 

9 

7 

325 

300 

327 

312 

8 

250 

259 

278 

276 

9 

110 

171 

186 

183 

10 

250 

259 

278 

276 

11 

58 

63 

77 

75 

12 

1194 

1280 

1313 

1300 

13 

148 

150 

162 

156 

14 

71 

76 

77 

77 

15 

21 

19 

23 

22 

16 

394 

415 

443 

423 

17 

235 

258 

265 

261 

18 

394 

415 

443 

423 

19 

272 

312 

328 

324 

cons_mths 

272 

312 

328 

324 

prec 

1 

91 

192 

204 

170 

2 

99 

206 

214 

248 

3 

104 

196 

202 

216 

4 

80 

192 

164 

172 

5 

93 

202 

188 

196 

6 

117 

266 

266 

256 

7 

148 

300 

324 

312 

8 

129 

264 

296 

278 

9 

98 

198 

204 

210 

10 

82 

152 

154 

154 

11 

71 

166 

172 

158 

12 

82 

226 

238 

230 

tmax 

1 

110 

104 

132 

107 
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cons_mths 


272 


300 


315 


342 
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Appendix  C:  Dataset  Generation  Processes 

To  ensure  that  they  had  complete  understanding  of  the  methodology  Cl  AT 
used  to  generate  their  data,  researchers  reviewed  the  available  on  line  lit¬ 
erature  and  summarized  the  processing  as  follows. 

C.l  WorldClim 

The  WorldClim  database  was  created  as  a  result  of  a  2005  study  published 
in  the  International  Journal  of  Climatology  (Hijmans  et  al.  2005).  Monthly 
averages  of  climate  at  weather  stations  were  compiled  for  a  period  of  ap¬ 
proximately  1950-2000  and  interpolated  using  a  thin-plate  smoothing 
spline  algorithm.  The  process  resulted  in  1  km  resolution  global  climate 
surfaces  (land  areas  only  -  excluding  Antarctica)  for  monthly  precipitation 
and  minimum,  mean,  and  maximum  temperature. 

Weather  station  data  were  collected  from  multiple  sources: 

•  Global  Historical  Climate  Network  Dataset  (GHCN) 

•  WMO  climatological  normals  (CLINO) 

•  FAOCLIM  2.0  global  climate  database 

•  International  Center  for  Tropical  Agriculture  (Cl AT)  database  assem¬ 
bled  by  Peter  G.  Jones  and  collaborators 

•  Additional  regional  databases  for  Latin  America,  the  Altiplano  in  Peru 
and  Bolivia,  and  the  “Nordic  Countries”  in  Europe,  Australia,  New  Zea¬ 
land,  and  Madagascar. 

As  a  measure  of  quality  control,  Hijmans  et  al.  used  ANUSPLIN-SPLINA 
4.3  to  implement  a  thin-spline  smoothing  procedure.  They  then  used 
manual  verification  of  spatial  and  elevation  data  for  weather  station  points 
in  conjunction  with  SPLINA  to  identify  the  stations  displaying  data  with 
the  largest  difference  between  the  “real”  data  and  that  in  the  spline  fitted 
climate  surface.  Data  records  thus  identified  with  errors  were  either  cor¬ 
rected  or  removed. 

The  construction  of  the  global  climate  surfaces  used  elevation  data  from 
both  the  Shuttle  Radar  Topography  Mission  (SRTM)  and  the  GTOPO30 
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database.  Monthly  surfaces  for  precipitation  and  minimum,  maximum, 
and  mean  temperature  were  created  while  uncertainty  in  the  interpolation 
process  was  dealt  with  through  two  strategies.  First,  the  authors  calculated 
elevation  bias  in  weather  stations  as  the  difference  between  the  mean  ele¬ 
vation  of  a  GCM  scale  grid  cell  and  the  mean  elevation  of  the  weather  sta¬ 
tions  therein.  Second,  they  used  SPLINA  to  build  a  surface  based  on  only 
half  of  the  stations,  then  mapped  the  differences  between  the  other  ob¬ 
served  half  and  the  predicted  values  to  illuminate  any  inconsistencies. 

To  illustrate  the  importance  of  high-resolution  datasets,  the  authors  calcu¬ 
lated  “within  cell  ranges”  for  a  10  arc-minute  resolution  grid  based  on  their 
l-km  resolution  surface  (10  arc-minutes  was  previously  the  highest  resolu¬ 
tion  global  climate  dataset  available  (New  et  al.  2002).  The  calculated 
ranges  of  annual  mean  temperature  and  annual  precipitation  were  com¬ 
pared  to  New  et  al.  (2002)  and  two  other  sets  of  high-resolution  climate 
surfaces  with  a  spatial  coverage  of  only  the  CONUS:  the  l-km  resolution 
Daymet  database  of  means  for  1980-1997  and  the  2.5  arc-minute  (-5km) 
PRISM  climate  database  for  1970-2000. 

Figure  Cl  shows  the  WorldClim  data  processing  flow  chart. 

C.2  International  Centre  for  Tropical  Agriculture  (CIAT) 

The  CIAT  spatially  downscaled  data  (Delta  Method)  (CCAFS  2011b)  is 
based  on  the  sum  of  interpolated  anomalies  to  the  high-resolution  month¬ 
ly  climate  surfaces  from  WorldClim  (Ramirez  and  Jarvis  2010).  The  da¬ 
taset  features  smoothed  (interpolated)  global  surfaces  of  multiple  (19-24 
depending  on  scenario)  GCMs  from  the  IPCC  Fourth  Assessment  Report. 

It  includes  GCM  output  from  three  emissions  scenarios  (SRES-AlB,  SRES- 
A2,  and  SRES-Bi)  for  several  30-year  running  mean  periods  (2010-2039 
[2025s],  2030-2059  [2045s],  2040-2069  [2055s],  2050-2079  [2065s], 
2060-2089  [2075s],  and  2070-2099  [2085s]). 


Figure  Cl.  WorldClim  data  processing  flow  chart. 
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Each  dataset  (comprised  of  an  SRES  scenario,  a  GCM  and  a  timeframe) 
includes  four  separate  monthly  variables:  mean,  maximum,  and  mini¬ 
mum  temperature,  and  total  precipitation.  Each  dataset  is  also  available  at 
four  spatial  resolutions  30  arc-seconds  (~i  km),  2.5  arc-minutes  (4.5  km), 

5  arc-minutes  (9.25  km)  and  10  arc-minutes  (18.5  km)  (1  minute  =  1.85 
km,  per  Penna  2002). 

The  Delta  downscaling  method  makes  the  following  gross  assumptions: 

•  Changes  in  climates  vary  only  over  large  distances  (i.e.,  as  large  as  the 
original  cell  size). 

•  Relationships  between  variables  in  the  baseline  (“current  climates”)  are 
likely  to  be  maintained  towards  the  future. 

Initial  data  were  compiled  from  the  Earth  System  Grid  (ESG)  online  plat¬ 
form  and  included  monthly  records  and  projections  from  1850-2100.  Fu¬ 
ture  projections  were  then  grouped  into  30-year  averages  before  anoma¬ 
lies  (delta  or  differences)  were  calculated  between  the  gathered  baseline 
climate  observations  and  future  projections  for  minimum/maximum  tem¬ 
perature  and  precipitation.  After  a  further  manipulation  involving  a  “thin- 
plate  interpolation”  in  which  the  anomalies  were  interpolated  between 
GCM  cell  centroids,  CIAT  was  left  with  a  30  arc-second  (1  km)  scale  sur¬ 
face  of  changes  in  climates  for  each  month  and  variable  in  the  dataset. 

The  final  interpolation  involved  combining  CIAT’s  downscaled  data  with 
that  from  WorldClim.  Interpolated  maps  were  added  to  the  “current” 
WorldClim  maps  (1950-2000)  account  for  the  possible  bias  due  to  the 
difference  between  WorldClim’s  downscaled  baseline  data  and  the  lower 
resolution  observations  obtained  through  ESG.  Finally,  mean  temperature 
was  calculated  as  the  average  between  minimum  and  maximum  tempera¬ 
ture  variables. 


Figure  C2  shows  the  CIAT  data  processing  flow  chart. 


Figure  C2.  CIAT  data  processing  flow  chart. 
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